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Isotropically Driven versus Outflow Driven Turbulence: 
Observational Consequences for Molecular Clouds 

Jonathan J. Carroll^, Adam Frank^, Eric G. Blackman^ 
ABSTRACT 



^' Feedback from protostellar outflows can influence the nature of turbulence in 

star forming regions even if they are not the primary source of velocity disper- 
sion for all scales of molecular clouds. For the rate and power expected in star 



forming regions, we previously ( ICarroU et al.l 120091 ) demonstrated that outflows 



Q> . coul d drive superson ic turbulence at levels consistent with the scaling relations 



C/^ ■ from iMatznerl (120071 ) although with a steeper velocity power spectrum than ex- 



r^ . pected for an isotropically driven supersonic turbulent cascade. Here we perform 

I . higher resolution simulations and combine simulations of outflow driven turbu- 

;_i . fence with those of isotropically forced turbulence. We find that the presence of 

<^ . outflows within an ambient isotropically driven turbulent environment produces 

a knee in the velocity power spectrum at the outflow scale and a steeper slope at 
sub-outflow scales than for a purely isotropically forced case. We also find that 
^ . the presence of outflows flattens the density spectrum at large scales effectively 

0> . reducing the formation of large scale turbulent density structures. These effects 

,__! . are qualitatively independent of resolution. We have also carried out Principal 

ly-^ . Component Analysis (PCA) for synthetic data from our simulations. We find 

O I that PCA as a tool for identifying the driving scale of turbulence has a mis- 

leading bias toward low amplitude large scale velocity structures even when they 
are not necessarily the dominant energy containing scales. This bias is absent 



K^ ! for isotropically forced turbulence but manifests strongly for collimated outflow 

driven turbulence. 

Subject headings: protostellar outflows, turbulence, star formation rate 

1. Introduction 

Supersonic velocity dispersions are a ubiquitous phenomena occuring over a wide range 
of scales in star forming molecular clouds. The relationship between velocity dispersions 
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(line-widths) and size was first sununarized by (lLarsonlll98ll ). Wh ile Larson's law s origi- 
nally applied only to global relationships between different CMC's, iHeyer fc BruntI (120041 ) 
showed that the same line-width size relationships held within individual CMC's. These ob- 
served velocity dispersions are commonly interpreted as resulting from turbulent gas motion. 
The formal distinction between turbulent motion and other sources of velo city dispersion in 



molecular clouds such as gravitational fragmentation is an important one (JField et al.l 12008 
20101 ). but for simplicity we presently use the terms turbulence and velocity dispersions 
interchangeably. 



The source of turbulent motions in molecular clouds is a topic of debate. If self-gravity is 
ignored then large scale sources of energy for turbulence include the kinetic energy injection 
from supernova and galactic differential rotation. Sources of turbulent energy on smaller 
scales include stellar feedback in the form of Hll regions, radiatively driven winds or accre- 
tion driven outflows. This raises the issue of the so-called injection scale, and whether there 
are in fact more than one dynamically important injection scale where energy is deposited 
into the turbulent spectrum. For a clump scale cloud (< lOpc), theoretical scenarios tend to 
divide the sources of turbulent forcing mentioned above into external and internal: External 
forcing models consider the injection scale to be larger than the cloud scale (e.g. Super- 
nova driven ISM turbulence) whereas inter nal forcing would refer to injection scales smaller 
than the cloud s ize (HII regions, outflows) (JMcKee &: Ostrikerl 120071 ) or winds from B stars 
(jValverde 1120101 ). Note that if self-gravity is included, then gravitational contraction from 
very large scales can itse lf produce velocity dispersions that are consistent with Larson's laws 
(JField et al.ll2008l . l2010l ). That represents an alternative paradigm for driving turbulent-like 
motions in molecular clouds. Whatever the source of the motions, in order for the turbu- 
lence in the parsec and subparsec scale clumps and cores to remain in an approximate steady 
state, it needs to be continually driven over time scales longer than a crossing time either 
internally via stellar feedback or externally via a cascade down from the ISM. 



Distinct from the question of the origin of turbulent motions, is the role of turbulent 'feed- 
back' on the star formation process. In the absence of turbulent support, most of the mass 
within a given structure (CMC, clump, core) woul d collapse into stars with in one free fall 
time giving an efficiency per free fall time e// = 1 (JMcKee fc Ostrikerl 120071 ). Observations 
however yield surprising low values of e// ranging from .01 to .1 not only for CMC's but 
also for substructures within CMC's of mean deii sities uh ranging from 10^ to 10^ cm~^ 
(JKrumholz &: McKed 120051 : iKrumholz fc Tan! 120071 ). One explanation for these low values 
of e// is the that the observed turbulent motions are signatures of internal kinetic support 
against collapse. We emphasize that the global energy budget of turbulence is a separate 
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issue from the the issue of local turbulent feedback. For example, star clusters form on scales 
that are a fraction of the total molecular cloud scale Ldoud/L duster ~ 10 and dynamical feed- 
back from outflows can impede star formation on cluster forming scales even if the total 
energy budget from outflows is insufficient to power the observed turbulent energy budget of 
the entire GMC This means that internal sources of turbulent forcing may be important on 
small scales whether or not they are important on large scales. On ev en smaller scales (ie. 
Lrnre ^ -1 Pc) outflows are likely responsib le for the dispersal of gas ( IVelusamy fc Langer 
19981 : iMatzner fc McKed 2000l: Myerd l2008l ) especially if the outflow opening angle widens 
with time ( Arce fc Sargentll2006l ). although the complex interaction between cores and their 
outflows requires more numerical modeling. In addition, not all forcings are equal: some 
forcings are isotropic and some are not. This distinction has observational consequences. 



In this paper we focus on the physics and observational consequences of coUimated outflow 
driving into both initially laminar and initially turbulent ambient environments. Theoret- 
ical treatments of feedback from mul tiple spherical o u tflows creating a self-regulating star 
forming systern was f irst e xplored by iNorman fc Silkl (119801 ) . Analytical work by Matzner 



flMatznei]l200d. 



2001 



20071 ) has explored the role of collimated outflow feedback on clouds. 



Krumholz et al.l (J2006l ) also considered the nature of feedback via outflows, concluding that 
these systems provide an important source of internal driving in d ense star f ormiri g cores. 
The flrs t simulation based stud y of the problem was prese nted by iMac Lowl (120001 ) . More 
recently iLi fc Nakamural (120061 ) and iNakamura &: Lil (120071 ) have mapped out the complex 
interplay between star formation outflows and turbulence. Their simulations include a self- 
consistent formulation of driving outflows from newly formed stars and they concluded that 
protostellar outflo ws were a viable means of generating turbulence in star forming clusters. 
Wang et al.l (120101 ) built on these studies using higher resolution AMR studies demonstrating 
the vigorous feedback that outflows are able to provide a cluster including the signiflcant 
modiflcation of the star formation efficiency and development of massive stars. We note 
that the Nakamura, Li and Wang studies all focused on cluster scales rather than molecular 
cloud scales. In addition they did not focus specifically on the mechanisms by which outflows 
generated tur bulence or the physic s of power spectra signatures. We also note a study of 
single jets by iBanerjee et al.l (120071 ) that explored the volume flUing averages of post shock 
concluding that single jets would not leave enough supersonic material in their wakes to act 
as a relevant source of internal forcing. 



The study of the detailed mechanics of outflo w driven turbulence was taken up in a series 
of paper s bvlCunningham et al.l (l2006al Jbl. l2008l ) . These investigations built off observational 
work by lQuillen et al.l ( 20051 ) for NGC 1333 which showed that turbulence was likely driven 
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not by active outflows bu t by the interaction of fossil shells driven by extinct outflows. In 
Cunningham et al.l ( 120081 ) it was shown how single outflows in a turbulen t environment would 
be disrupted leading to their energy resupplying turbulent motions. In ICarroU et al.l (|2009[ ) 
a study of multiple, randomly oriented jets driven into a quiescent medium was undertaken. 
The results of these simulations showed conclusively that interacting outflows can lead to 
turbulent motions. Explorations of the energy spectrum and the density PDF showed that 
statistically steady isotropic random motions were maintained by the outflows. Of particular 
interest was the fact that an injection scale associated with outflow cavity interaction Koutfiow 
was apparent in the energy spectrum E{k). Also the power law index of the spectrum was 
steeper than would be expected for Burgers type turbulence su c h as t hat generated when less 
reahstic isotropic forcing at large scales is used. ICarroU et al.l (120091 ) found that the steeper 
spectrum was associated with outflow cavities sweeping up vortices at smaller scales. Thus 
the Carroll results provided hath proof that outflows can drive turbulence and provided some 
insights into the physics involved and its statistical signatures. 



The goal of discerning scales of tu rbulent driving in real molecular clouds has, recently, 
received considerable attention. In JBrunt et al.l (12009! ) PCA was used to determine the 
dominant scales of turbulent driving based on molecular line emission. Based on comparisons 
of first and second eigenimages of molecular channel maps for ^^CO and ^^CO they find that 
the data support the conclusion that that the bul k of the energy driv ing for turbulence 
comes from scales much larger than that of clumps. JPadoan et al.l (120091) re ached a similar 
conclusion by appl ying VCS (Velocity Coordinate Spectrum) (lLazarianll2009l ) on ^'^CO data. 
Brunt et al.l (120091 ) also analyzed C^'^0 data which probes denser regions of the clumps, and 
this data did reveal evidence for a second, subparsec driving scale consistent with that of 
outflows. While neither of these studies preclude the importance of feedback from outflows 
on very small scales, neither of these papers actually invoke full 3-D simulations of coUimated 
outflow driven turbulence in their conclusions about driving scale. Thus there is a potential 
loophole in the application of methods that require a density-velocity scaling relation in 
order to interpret observations as a signature of kinetic energy density on a given scale. One 
of the topics we address in this paper is the bias that results from this specific issue. 



In section 2 we discuss the basic parameters and setup of our simulations and choices of 
forcing for both our coUimated outflow driving and our isotropically forced cases. In section 3 
we discuss the results comparing the isotropically driven turbulence with coUimated outflow 
driven turbulence and the interplay between the two. In particular, we analyze the growth 
and saturation of turbulence, the velocity and density spectra, the characteristic driving 
scales, and the relative insensitivity to resolution. In section 4 we uncover biases in the 
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current use of PCA to infer the scales of turbulent driving when comparing isotropically 
forced to outflow driven turbulence. We show that PCA fails to identify where the dominant 
driving power is when comparing sources of turbulence with different forcings. We conclude 
in section 5. 



Numerical Model 



As in lCarroU et al.l ( 120091 ) we used an MHD code called AstroCUB but without magnetic 



fields for this work. The code is 2nd order accura te in space and time and uses a non-split 



CTU method with upwinded CT as described by iGardiner fc Stond (120081 ). All four of the 
simulations were performed on a periodic cube of length Lbox = 1-48 pc initialized with a 
uniform density p^ = 2.51 x 10~^''g cm~^ and mean particle weight of 2.1 amu at 10 K giving a 
sound speed of .20 km/s. Each of the simulations was driven by some combination of outflow 
feedback ( §2.ip and/or isotropic forcing ( §2.2p described below. Run HDO was forced solely 
via outflows, run HDI was forced isotropically, and run HDOI employed both forms of forcing 
each of comparable strength. It is worth mentioning that the difference in velocity between 
typical outflows (240 km/s) and the turbulence they generate (1 km/s) makes simulations 
that model both types of flows computationally intensive. This is because with Eulerian 
grid-based codes, the computational time steps are hmited by the highest speeds (240 km/s) 
but the timescale of interest is set by the turbulent velocity (1 km/s). Therefore since 
simulations with outflow forcing would normally take of order 100 times as long as those 
with just isotropic forcing, we started runs with outflow forcing at lower resolutions and then 
periodically doubled the resolution until the desired resolution was reached. Since supersonic 
turbulence decays on the order of 1 dynamical time t^ = -^^^^, any effects from regridding are 
expected to also disappear after ltd- For this reason, all of the spectra were taken after at 
least ltd after the final regridding. All of the runs and the regridding times are summarized 
in Table [H 



2.1. Outflow Feedback 

During the runs with outflow feedback, highly collimated bipolar outflows of momentum 
V = 20.0Mq km/s were launched at the volumetric rate oi S = 58.4pc^'^Myr^^ or with a 
period of Tiaunch = SLil^^ = 5.28 kyr. Combining the outflow momentum V and rate S with 
the mean dens ity p^ determines the characteristic outflow scales of mass, length, and time 
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Combining these gives other characteristic quantities. Of particular interest is the charac- 
teristic velocity, the characteristic wave number, and the characteristic accleration: 

Table [2] summarizes the various resulting outflow scales. 



Each outflow was launched from a cylindrical source region randomly located and oriented 
with radius r^ and length 2zo. Instead of constantly setting the values of density and mo- 



mentum within each launch region as in ICarroU et al.l ( l2009l ) , source terms were calculated 



to supply the mass and momentum needed to maintain the desired density and velocity pro- 
files. The density profile was chosen to be constant {p^^) while the velocity profile contained 
a quiescent core {\z\ < Zi) with \z\ the direction of the outflow. Beginning with bipolar 
acceleration regions from Zi < \z\ < Zo the proflle is described below. 

_ ■.Q<\z\<Zi 

' sign {z) °^'^'^ '-' ■■ Zi<\z\<Zo 

The source terms for density Sp and momentum Sp^ were calculated by substituting the 
launch proflles into the conservation equations and solving for a steady state solution. 

|e = _v ■ (pv) + 5, = ^ 5p = V ■ (pv) (1) 

^ = -V ■ (v : (pv)) - VP + 5,v = ^ 5,v = V ■ (v: (pv)) + VP (2) 

As was mentioned before, modehng outflows and the turbulence they produce is computa- 
tionally difficult due to the high ratio between the outflow velocity V^ and the turbulent 



Table 1. Table of Runs 

Run Isotropic Outflow Initial Final Regrid Time in T 

Forcing Forcing Resolution Resolution 128 256 512 

5 6 



HDO 


No 


Yes 


128^ 


5123 


HDI 


Yes 


No 


512^ 


5123 


HDOI 


Yes 


Yes 


64^ 


5123 
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velocities they produce V. Fortunately, since outflows are highly radiative, it is their mo- 
mentum which is most critical in modeling their feedback. This allows us to increase the 
mass injected while decreasing the velocity thereby reducing the ratio V^/V. So while the 
momentum injected by each outflow was modeled on a .SM© star ejecting 1/6 of its mass at 
240 km/s, the actual mass injected by each outflow was increased from .OSSMq to .305Mq 
to allow for a slower velocity jet (V^/V = 61 instead of 224). This outflow mass loss increases 
the mean density by ~ 2% for each outflow time T. Since the simulations were run for 8T 
the overall increase in the mean density was < 16%. In addition the mass loss rate was 
chosen to be 1.30 x IO^^Mq/jt giving an outflow duration of 2.35 kyr. The actual duration 
of the outflow was modified slightly to allow for a gradual ramp down of the velocity while 
keeping the total momentum injected equal to V. This helped to prevent a large rarefaction 
from developing once the outflow sourced terms were shut off. Since the outflow period was 
Tiaunch = 5.28 kyr, only one outflow was ever active at a given time. Table [3] summarizes the 
outflow parameters. 



2.2. Isotropic Forcing 

In addition to forcing via outflow feedback, a constant isotropic solenoidal forcing at the 
box scale L^ox = 4£ was used to mimic the cascade of turbulent energy from larger scales. 
The mean acceleration of the isotropic forcing was (a) = .95^ while the components and 
phases of each wavevector were chosen randomly. Table H] lists the forcing vectors and phases 
for each wave mode. The total force applied to each cell was 

f = pa{x)=p^AnCOs(kn-X + (l)n] (3) 

n 

In a similar manner to the outflow forcing, we can combine the isotropic forcing scale 
Lbox with the mean acceleration A and the ambient density Pg to calculate characteristic 
scales for the isotropically driven turbulence. Of particular interest is the velocity scale 
V/ = ^yLboxA. = ViJCA = 2V and the associated time scale 7/ = Vi/A = 2T. While both 
types of forcing have the same mean accleration and inject momentum at comparable rates, 
the isotropic forcing has a greater velocity scale simply because it injects energy on a larger 
spatial scale. For completeness, all of the equivalent scales associated with the isotropic 
forcing are listed in table [2] 



Table 2. Outflow Scales 





cgs units 


astronomical units 


isotropic scale 


Po 


2.51 X 10-20 g cm-3 


371 M© pc-3 


Po 


V 


3.98 X 10^9 g cm s^^ 


20.0 Mq km/s 


128P 


s 


6.31 X 10-6S cm-3 s-i 


58.4pc-3 Myr-i 


5/128 


M 


3.73 X 10^^ g 


18.7 Mo 


UM 


C 


1.14 X 10^^ cm 


.370 pc 


AC 


r 


1.07 X 10^3 g 


.338 Myr 


2r 


V 


1.07 X 10^ cm/s 


1.07 km/s 


2V 


A 


1.00 X 10-^ cm/s^ 


3.23 pc/Myr2 


A 



Table 3. Outflow Params 



description 


symbol 


value 




comment 


density- 


Po 


2.5 X 10-20 g/cm^ 




velocity 


Vo 


65.5 km/s 






duration 


to 


2.34 kyr 






period 


To 


5.28 kyr 






radius 


To 


5960 AU 




lOAx at 512^ 


inner buffer 


Zi 


2Ax 




1190 AU at 512^ 


launch thickness 


Az 


3580 AU 




6Aa; at 512^ 


opening angle 


Oo 


0° 






mass loss rate 


M 


1.30 X IO-^Mq 


yr-i 
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3. Results 

3.1. Growth and Saturation of Turbulence 

Figure [T] shows the growth and saturation of the mean scalar momentum P =< p\v\ > 
and the velocity dispersion frms = v<v^> with just isotropic forcing (HDl), just outflow 
forcing (HDO), and where both types of forcing are present (HDOl). Since both types of 
forcing have approximately the same mean acceleration [A for the outflow forcing, and .95^ 
for the isotropic forcing), the initial growth of the scalar momentum seen in the left panel 
of flgure [1] for run HDOl is larger (about twice) that of runs HDO and HDl as expected, 
although due to the lower initial resolutions of run HDO and HDOl, the initial growth in the 
outflow momentum injection rate appears to be ~ 80%PgV/T. The effects of the resolution 
can also be clearly seen in the resaturation of the scalar momentum for run HDO following 
regridding at t = 5T and 6T and somewhat less so for run HDOl at t = 3T, 5T, and 6T. 
Each time the resolution is doubled, the dissipation scale is halved, the dissipation rate is 
reduced, while the injection rate remains the same - leading to a higher saturation level. 

Even though both types of forcing have the same mean acceleration, the flnal saturation 
levels for the mean scalar momentum for runs HDl and HDO differ due to the difference in 
length scales associated with the driving. The isotropic forcing injects its energy/momentum 
at a larger scale {Lbox = 4£ instead of £), leading to a longer turnover or cascade time (2T 
instead of T), and a slower dissipation. This results in a higher saturation level consistent 
with the velocity scale for the isotropic forcing (V/ = 2V). The saturation level for run 
HDOl is slightly larger then that of run HDl due to the additional momentum injected by 
the outflows. Given the scaling relations, we might expect the velocity dispersion for runs 
HDl and HDO to follow a similar trend. Instead, we see that in spite of having a lower 
turbulent velocity scale, the outflow driven turbulence has a comparable velocity dispersion. 

Table 4. Isotropic forcing components 



k A 



[1,0,0] [0, .04087, .81547] 1.6519 

[0,1,0] [-.64450,0, -.50128] 1.1919 
[0,0,1] [-.41577, .70271,0] 3.7784 



-10- 



This can be explained by considering the high velocity nature of the outflows themselves 
(Vo/V = 61). It takes very little high velocity gas to significantly contribute to the velocity 
dispersion, and when outflows are present there is plenty of it around. In addition, the high 
velocity impulsive nature of the outflow driving leads to the high degree of variability in the 
velocity dispersion. 

In flgure [2] we show joint probablility distribution for the three hydodynamic runs. The 
isotropically forc ed runs show' s the e xpected log- normal distribution with density. As was 



demonstrated in ICarroll et al.l (120091 ) outflows will also produce such a log-normal distribu- 
tion in p as can be seen in the higher resolution HDO runs. Thus we conflrm the conclusions 
that transient outflow cavities can set the bulk of initially quiescent material into random 
but statistically steady supersonic motions. The pdf's for runs HDO and HDOI also re- 
veal the presence of the high velocity material produced by the outflows responsible for the 
high velocity dispersions and the higher degree of variability. The isotropic forcing on the 
other hand is only able to accelerate material at A over a distance Li,ox to velocities of 
~ yj2Lbox'A. = \/8V leading to the cutoff in the left panel of flgure [2] (run HDO) and the 
horizontal kink in the middle panel (run HDOI). 



3.2. Driving Scales of Turbulence 

In order to better understand how these two types of forcing influence their environments 
over a range of scales, we looked at the velocity and density spectra for runs HDO, HDI, 
and HDOI after a statistically steady state had been reached at a resolution of 512^. First, 
as would be expected, the driving scales associated with the outflow forcing (/C) and the 
isotropic forcing (-r^— = 1) are clearly evident in the velocity spectra seen in flgure O 
The velocity spectra for runs HDO and HDOI have a clear peak and knee respectively 
at the outflow wave number /C, while the velocity spectra for runs HDI and HDOI peak 
at the isotropic forcing scale -t^- = 1. Secondly as would also be expected, run HDI has a 
spectral index of /3 = 2.0 over the marginally resolved inertial range, consistent with previous 
simulations of supersonic, isothermal, isotropically-forced turbulence as well as supersonic 
cascade and shock-dominated models. Somewhat unexpectedly, however, is the steep slope 
{(3 = 3.2) present in the velocity spectra for runs HDO and HDOI at sub-outflow scales. 
This steep slope along with the signiflcantly greater amount of energy imply more coherent 
higher-velocity structures on sub-outflow scales. This is reasonable however, given the nature 
of the outflow forcing. Outflow cavities are able to expand coherently out to distances £ 
before they have decelerated to the background turbulent velocities, at which point they 
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Comparison of Mean Scalar Momentum 



Comparison of Velocity Dispersion 




4 6 

Time {T) 




Fig. 1. — Time series showing development of the mean scalar momentum density P and 
the velocity dispersion Vmis for runs HDI (solid lines), HDOI (dash-dotted lines), and HDO 
(dashed lines). All quantities are scaled to the outflow scales (ie T,V, & Po)- 



p-v pdf for run HDI 



p-v pdf for run HDOI 



p-v pdf for run HDO 



2 

1.5 

1 

0.5 

^ -0.5 
-1 

-1.5 
-2 

-2.5 



Contour levels are 
separated by factors of 2 





-2 

log(p/A,) 




-4-2 2 -6 

log(p/A,) 



-4-2 2 

log(p/A,) 



Fig. 2. — Joint probability distribution functions of velocity (y-axis) and density (x-axis) 
for runs HDI (left), HDOI (middle), and HDO (right). Lighter areas correspond to higher 
probability. The crosshairs are at the mean density Pg and the outflow scale V 
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lose coherence. As the cavities break up they will stir motions with velocities V over size 
scales L and inject their remaining momentum into turbulent motions. The time scale for 
these turbulent motions to cascade however y = T, is the same as the timescale between 
successive outflows, so each outflow effectively destroys any decoherent motions left by its 
predecessors. It is interesting that the source of the turbulent energy is also, to some degree, 
the source of its dissipation. We conjecture that this conclusion will hold for all forms of 
internal driving where stellar sources feed energy back into the parent cloud. Our results 
highlight the difference between simulations which track the more realistic discrete driving 
and those which impose an apriori driving spectrum at large scales to model an unknown 
process of external driving. 



In contrast to the effect outflows have on the velocity spectra (a steepening at sub-outflow 
scales), the density spectra undergoes a flattenning at supra-outflow scales as seen in the 
right panel of figure El Note the density spectrum for run HDO is quite flat at large scales 
and that run HDOI is flattened compared to run HDL This is not too surprising when one 
considers that each outflow is randomly placed and oriented. To gain a visual sense of what 
this means we present in the top panels of figure [5]2D cuts of density through both the HDOI 
and HDI simulations once saturation has been achieved. In addition, in the bottom panels 
of figure we show log column density maps to give a sense of the difference in observational 
appearance the two modes of turbulent driving would obtain. 

When only outflows are present (as in run HDO), there is no source for large (supra- 
outflow) scale coherent motions producing large scale coherent density structures. The ran- 
dom placement and orientation of outflows produces fairly random density structures and 
a fairly white spectra. When only large scale isotropic forcing is present (as in run HDI) 
do large scale motions driven by large scale isotropic forcing generate large scale coherent 
density structures. These can be clearly seen in the right panels of figures [51 The time scale 
however needed for these structures to form is 7} = IT . Thus when both large scale forcing 
and outflows are present (as in run HDOI), these large scale density structrues are unable 
to survive being torn apart by sequential rounds of outflow generation (which occur on a 
timescale T). Thus outflows suppress small k modes in density resulting in a flatter density 
spectra as seen in figure [31 as well as a less coherent density field as seen in the left panels 
of figure [51 This disruption in large scale density structures is somewhat mirrored in the 
disruption of large scale velocity structures. This explains the decrease in power at large 
scales seen in run HDOI as compared to run HDI. 
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3.3. Resolution study 

To better understand the role numerical dissipation might play in the density and ve- 
locity spectra, we completed run HDOl at 64^, 128^, 256'^ as well as 512^ in order to perform 
a resolution study. Figure H] shows both the velocity and density power spectra at ~ 8T 
at resolutions of 64^, 128^, 256^, and 512^. The velocity spectra are fairly insensitive to 
numerical resolution, especially at sub-outflow scales. In cascade models of turbulence, the 
lack of turbulent eddies near the dissipation scale reduces the rate at which energy can cas- 
cade creating a pileup of energy at scales above the dissipation scale. This bottleneck and 
subsequent steepening depends entirely on the dissipation scale and would be resolution de- 
pendent. The lack of any resolution dependent bottleneck effect is however consistent with 
the lack of any cascade at sub outflow scales, owing to outflows themselves being responsible 
for the turbulent dissipation. The outflows themselves create a kind of bottleneck, but this 
bottleneck is at the outflow scale - not the dissipation scale - and is well resolved (except 
perhaps marginally so at 64'^). 

The density spectra, on the other hand, due show strong resolution dependent effects. At 
supra-outflow scales, the spectra become flatter as the resolution increases due to the ability 
of the outflows to more flnely shred the large scale density structures. The spectra also 
appear to steepen at a resolution dependent scale of ~ 8Ax most likely due to numerical 
diffusion near shocks. At a resolution of 512^, this steepening scale is a few times smaller than 
the sonic length l^ ~ 24Ax where v{ls) = Cg or P^,(/s) = c^ = .035V^. Below the sonic length 
scale, turbulent velocity fluctuations become subsonic (and therefore less compressible). 



3.4. Discussion of Dynamical Results 

Outflow driven turbulence is very different from isotropically forced turbulence, and the 
presence of outflows within an istropically driven cascade signiflcantly alters the dynamics. 
In general the presence of outflows will lead to a steeper velocity spectrum at sub-outflow 
scales and a shallower density spectrum at supra-outflow scales. The steepening of the 
velocity spectrum at sub-outflow scales can be attributed to the presence of coherent outflows 
themselves, as well as their ability to effectively impede the cascade of energy by sweeping up 
small eddies as they expand. As a result, the dynamics on sub-outflow scales are dominated 
by the passing of subsequent outflows rather than the cascading of energy from larger scales. 
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The flatter density spectra at supra-outflow scales can be attributed to the ability of the 
outflows to disrupt larger scale structures from forming. This can be most clearly seen upon 
visual inspection of the bottom panels in figure O When the outflows are highly coUimated 
(as was the case here), they can also disrupt the coherent density structures produced by 
neighboring outflows leading to a flattened density spectra even into the sub-outflow scale as 
seen in figure |3l Outflow feedback is therefore important in not only sustaining turbulence, 
but also in limiting the spatial scales of density structures and therefore shaping the IMF 
(initial mass function). This interplay will be somewhat modified when the outflow locations 
and orientations are linked to the glo bal dynamics of th e clouds via gravitational collapse 



and magnetic fields, but we see as in IWang et al.l (120101 ) that outflows have the potential 



to play a large role in setting the star formation rate and the initial mass function in their 
neighborhoods. 

We note that our results do not imply that outflows can drive turbulence across the entire 
cloud. Given the paucity of YSO outflows across the entire length of a cloud like Perseus 
it seems unlikely that the cloud as a whole can be set into turbulent motions from outflow 
feedback. But this does not preclude outflow driving in regions where stars actually form. 
The pervasive clustering inherent to star formation implies that outflows can significantly 
affect the star formation rate even though they likely do not drive the global GMC turbulence. 



4. Connecting to Observations 

Although the model of outflow driven turbulence used here is somewhat idealized by 
the outflows having uniform strength, collimation, and mean separation, it is reasonable 
to expect a similar signal to be present in the actual velocity spectra for regions where 
outflows are dynamically important. We note that such a signal is likely to be blurred since 
outflows will have varying strength, collimation, and separation. Even under the best of 
conditions, observers are still left with the difficult task of inferring the velocity spectra 
using VGA (velocity channal analysis), VGS, or PGA from data cubes of emission intensity 
as a function of position on the sky and line of sight velocity T{x, y, Vz)- Since emission is at 
least density dependent, some assumptions on density- velocity correlations and distributions 
must be made. These assumptions are often justified on the basis of both theoretical and 
numerical models of isotropically forced turbulence. Figure |2]however, shows clearly that the 
density-velocity distributions and to some extent correlations are different between outflow 
driven turbulence and isotropically forced turbulence. It is therefore worth investigating the 
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ability of various techniques to reconstruct the velocity spectrum of models of outflow driven 
turbulence (as well as more realistic simulations that involve self-gravity etc.). Here we look 
at the efficacy of one such technique (PCA) to infer the presence of outflows. 



4.1. Synthetic Observations 

In practice, the actual emission from molecular clouds will depend on more than just 
density (ie temperature and chemistry). In addition, optical depth effects further complicate 
the issue. Here we ignore such effects and assume ideal observations of a density weighted 
velocity field. This allows us to focus more directly on how the different density-velocity 
correlations and distributions produced by isotropic and outflow forcing effect the results 
inferred from PCA. To convert our computational box into a spherical cloud we first tiled 



the bo x in each direction. This allowed us to apply a spherical window as in iBrunt et al. 



(120091 ) without discarding as much information. This also increased the ratio of cloud size to 
isotropic forcing scale and prevented a sort of blending of the first two principal components 
(described below). We then created synthetic data cubes created by making density weighted 
histograms of the line of sight velocity along the x-axis at each yz position for runs HDI, 
HDOI, and HDO. 



4.2. Global Line Profiles 

Before applying PCA it is instructive to examine the 'unresolved' spectra (integrated 
over the plane of the sky) for each run (upper left panel of figure [6]). Note that runs with 
isotropic forcing (HDI and HDOI) have a broader central line profile then the run with just 
outfiow forcing (HDO) roughly consistent with the greater velocity scale V/ = 2V associated 
with the isotropic forcing as well as the saturation levels of momentum seen in the left panel 
of figure [H It is not surprising that the linewidths are more consistent with the saturation 
levels of momentum as opposed to velocity dispersions since both momentum and emission 
are weighted by density. 

Additionally the line profiles for runs with outfiow forcing (HDO and HDOI) have extended 
wings not present in run HDI due to the higher velocity material injected by the outfiows. 
While the wings are only marginally apparent in the density weighted emission and as such do 
not contribute to the observed line- widths, they do contribute significantly to the momentum 
and energy. This becomes apparent after weighting the 'observed' line profiles by the line 
of sight velocity v^ or by f^. These modified spectra seen in the bottom panels of figure 
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|6] show the distribution of momentuin or energy at the different velocities. Note the large 
contribution to the momentum and energy provided by the high velocity wings in runs HDO 
and HDOI. As a result, line-widths inferred from density weighted emission will tend to 
underestimate the energy and momentum injected by the outflows. In addition to the line 
profiles produced by density weighted emission, we constructed line profiles of volume (or 
unweighted) emission (upper right panel of figure [6]). Note the global line profiles are very 
similar implying a lack of any strong density velocity correlations. This is not surprising 
since the density- velocity distributions in figure H] are only weakly correlated over the bulk 
of the material for all three runs. 



4.3. Principal Component Analysis 

Principal Component Analysis is a technique for finding relationships between a given 
set of variables that can best explain the variance in a given set of measurements. Given a 
set of n measurements of p variables Tjj (corresponding to the i^^ measurement of variable 

n 

j), PCA first involves the construction of the covariance matrix Sjk = 2_,TijTik where we 

n 

have subtracted off the mean values for each variable (Tjj = Tjj — Tj where Tj = ^ y ^Tjj). 



i=l 



The first principal component is then the (normalized) eigenvector V^ of the covariance 

p 

matrix (and corresponding projection P/ = /^^jj^-^) with the largest eigenvalue Aj and as 

i=i 
such can account for the largest source of variance. The second principal component is the 

eigenvector with the second largest eigenvalue and so on. 



Principal Component Anal ysis as applied to interst ellar turbulence is formally described in 



Heyer fc Schloerbl (119971 ) and iBrunt &: Heyerl (l2002al ) . It involves treating the spectroscopic 



datacube Tjj = T{xi,yi,Vj) as a measurement at each pixel {xi,yi) of a set of variables (the 
emission in each velocity channel Vj). The eigenvectors of the covariance matrix are then 
line spectra in velocity space, and the projections are images in xy space where each pixel 
represents the degree of correlation between the line spectra at that pixel and the principal 
eigenvector. These principal 'eigenspectra' and corresponding 'eigenimages' then account 
for the largest sources of variance in the datacube. While PCA formally uses the covariance 
matrix, a modified version (mPCA) has also been used to decompose datacubes into eigen- 

n 

spectra and eigenimages - though using the 'co-emission' matrix Sjk = y^TijTik in which 
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the mean values for each channel are not subtracted. The ij^'^ element of the co-emission 
matrix is therefore the degree of correlation between the emission in velocity channels i 
and j summed over all positions. The lack of any mean value subtraction introduces addi- 
tional pseudo-variance which is effectively remove d in the first principal component. Here 



we use mPCA to make contact with previous work (JBrunt fc Heyeiil2002al : iBrunt et al.ll2009 



Brunt fc Hevei]l2002a[ l 



PCA effectively reduces the dimensionahty of the dataset by projecting the 3-D datacube 
onto a sequence of 1-D principal eigenspectra resulting in a sequence of 2-D principal eigen- 
images. The generation of the eigenspectra (and therefore the eigenimages) however, does 
not make use of the velocities associated with each channel (or the positions on the sky of 
each pixel). The positions in the sky of each spectra could be shuffled without changing 
the eigenvectors, and the velocity channels could be shuffled without changing the resulting 
eigenimages. However, the resulting (unshuffled) eigenvectors and eigenimages can be used 
to determine sequentially the velocity magnitudes {dvi) and length scales (/j) over which 
the largest sources of variance change. These scales are typically calculated as the distance 
at which the autocorrelation functions of the eigenspectra (dvi) and eigenimages (/«) drop 
by a factor of e. The combination of /, and dvi for each principal component then allow 
for the reconstruction of a line-width size relationship (or equivalently a velocity spectra 
^pca)- This inferred spectra has been sh own to closely rnimic t he actual veloci ty spectra in 



fract ional Brownian motion simulations (JBrunt &: Heyerll2002al ). In addition, (JBrunt et al. 



20091 ) demonstrated that the ratio of I2 to h (when using mPCA) is correlated with the ratio 
of cloud size Xd to driving scale L^ox in simulations of isotropically forced turbulence. Since 
the first principal component in mPCA is likely to mirror the mean values to account for 
the pseudo- variance, the principal length /i is typically the size of the cloud. This leaves the 
second principal component and corresponding length I2 to account for the largest source of 
true variance which should be associated with the driving scale. 



Figures [7] and [8] show the first three eigenimages and eigenvectors for the three different 
cases. Note that while all three runs have similar velocity dispersions (figured]), the runs 
with isotropic forcing (HDI and HDOI) have much broader eigenspectra (figure [H]). This 
broadening results in higher spectral correlation magnitudes dvi for runs HDI and HDOI 
given in table [51 These differences however, are in general consistent with the differences 
seen in the final values of the mean scalar momentum in figure [H and is not surprising 
considering that both the mean scalar momentum and emission are density weighted. It is 
also interesting to note that the first three spectral correlation lengths for run HDOI are 
actually smaller {^ 10%) than those for run HDI even though run HDOI has additional 
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forcing via outflows, a larger velocity dispersion, and a larger saturation value of mean 
scalar momentum. We interpret this decrease in velocity correlation length as arising from 
the disruption by outflows of the large scale coherent flows produced by the isotropic forcing. 
This is consistent with the decrease in power at large scales seen in the velocity spectra in 
figure [31 

While the presence of outflows in run HDOI does not contribute to larger velocity corre- 
lation lengths {dvi, dv2, & dv^) and as such is not readily apparent in the principal eigen- 
spectrum, the presence of outflows is apparent in the first principal eigenimage in the lack of 
coherent large scale density structures for run HDOI compared to run HDL This is consistent 
with the flattening of the density spectra seen in figure [3] and the larger principal autocorre- 
lation length (/i) given in table [5l However, the second and third eigenimages for run HDOI 
are very similar to run HDI and show clear evidence for the large scale driving. These large 
dipole structures present in the second eigenimages for runs HDI and HDOI result in larger 
autocorrelation lengths {I2) a nd larger ratios of /2 to /i especially when compared to run HDO 



see table [5]). We find as in (JBrunt et al.l 120091 ) that the ratio of I2 to h correctly identifies 



the largest driving scale. Runs HDI and HDOI are driven isotropically with y^- = 1 and 
have ratios f- = .3649 and .2723 respectively, while run HDO is driven by outflows at one 
fourth the scale {j^- = -25) and has a corresponding ratio roughly one fourth (^^ = .0874. 
While the ratio of I2 to li does not reveal the presence of the outflows in run HDOI, the 
resulting line- width size relationships shown in the left panel of figure [9] do show an excess of 
energy for run HDOI beginning at the third principal component whose correlation length 
^3 = .1142 is approximately the same as that for the second principal component for run 
HDO I2 = .0973 which is the scale associated with the outflow driving. We note however, 
that since the inferred spectra from PGA relies on line widths that are fairly insensitive to 
the high velocity wings, PGA will in general underestimate the true strength of the outflow 
forcing. This explains why run HDI appears to have more power over the entire range of 
scales in the results from PGA when in fact run HDI has less power at small scales as seen 
in the true velocity spectrum in figure [31 



We have already seen how the density weighting of the emission results in linewidths 
consistent with values of the mean scalar momentum rather than the velocity dispersion. 
This alone explains why the second principal component for run HDOI picks out the variance 
produced by the isotropic forcing rather then that due to the onflows. However, since the 
emission is not weighted by velocity, low velocity motions on a given scale can produce just as 
much variance in the datacube as high velocity signals on the same scale. Thus a low velocity 
large scale coherent motion could produce more variance than a high velocity motion of less 



-19- 



coherence on that scale. Thus the ratio of I2 to li is more a simple measure of the largest 
scale of coherent motion rather than the scale of dominant driving energy. In simulations 
of isotropically forced turbulence there is no distinction. But in clouds where outflows are 
present within an external cascade, or where there are large scale coherent motions due to 
angular momentum conserving collapse, this distinction becomes important. 

To further test the aforementioned bias we ran PCA on a series of datacubes produced 
by modifying the velocity field in run HDO through the addition of a large scale sinusoidal 
variation in the line of sight velocity - similar to what is seen in run HDOI. We varied 
the strength of the perturbation from 10%V to 50%V in 10% increments (HDOIO, HDO20, 
HDO30, HDO40, & HDO50) and looked at how the resulting correlation lengths shifted. 
The first three eigenimages for each run are shown in figure [TOl By run HDO30, the second 
principal eigenimage shows the large scale velocity variation and the resulting autocorrelation 
length I2 has increased from its initial value of .0587 to .1477 corresponding to the driving 
scale seen in runs HDI and HDOI. Note also the (small) 20% increase in dvi from .10245 
to 1.2009 produced by the velocity perturbation. These trends are also evident in the line- 
width size relationships shown in the right panel of figure HI Thus our results indicate that 
PCA is not a reliable measure for determining the specific question of the dominant mode 
of energy input in driving turbulence. In our experiments even though the energy put into 
the solenoidal motions was smaller than that in the outfiows, PCA still picked up these 
solenoidal motions as dominant i.e. Esoi < Eof but l2{soi) > h{of)- 



5. Conclusion 



Building on the work of ICarroU et al.l (120091 ) we studied the properties of turbulence 



driven isotropically, by outfiows, and by a comb ination of bot h for in put parameters typical 



of clump scale molecular clouds. We find as in ICarroU et al.l (|2009[ ) that outfiows are able 



to drive supers onic turbulence at levels consistent with the scaling relations put forth by 



Matznerl (120071 ). We also find that turbulence driven by outfiows is characterized by a steeper 
velocity spectrum at sub-outfiow scales, and an overall fiatter density spectrum. Based on 
our resolution study, we conclude that when outfiows are present, the resulting steep velocity 
spectra and correspondingly smooth velocity structures seen at sub-outfiow scales are not 
due to numerical dissipation, but rather due to the presence of the outfiows themselves. 
Outfiows are able to continually generate smooth velocity structures on outfiow scales that 
sweep up and 'dissipate' smaller structures at the same rate that these structures would 
begin to cascade. Outfiow driven turbulence is therefore quite different from the cascade of 
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Fig. 3. — Velocity power spectrum (left) and density spectrum (right) for runs HDI (solid 
line), HDOI (dash-dotted line), and HDO (dashed line). The vertical dashed line corresponds 
to the outflow wave number /C and all quantities are scaled to the outflow scales. 



Table 5. Autocorrelation lengths of eigenimages and eigenvectors from principal 
component analysis. All lengths are in pc and velocities in km/s 



Run 


k 


k 


k 


kik 


^o/Lhox 


dvi 


dv2 


"^rms 


f^rms 


HDI 


.7960 


.2904 


.2076 


.3649 


1.000 


2.376 


.9276 


1.9785 


1.201 


HDOI 


1.126 


.3067 


.1142 


.2723 


1.000 


2.183 


.8593 


2.2546 


.9682 


HDO 


1.113 


.0973 


.0492 


.0874 


.2500 


.7977 


.3293 


2.2907 


.3482 


HDOIO 


1.115 


.1251 


.0537 


.1123 




.8187 


.3289 






HDO20 


1.117 


.1964 


.0595 


.1758 




.8684 


.3520 






HDO30 


1.119 


.3138 


.0703 


.2804 




.9402 


.3903 






HDO40 


1.118 


.3874 


.0927 


.3466 




1.023 


.4370 






HDO50 


1.112 


.4367 


.1235 


.3928 




1.108 


.4907 
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Velocity Spectra Resolution Comparison 




Density Spectra Resolution Comparison 




Fig. 4. — Plots showing velocity and density power spectrum for run HDOI at different 
resolutions. Note the insensitivity of the velocity spectrum to resolution at sub-outflow 
scales. Also note the strong resolution-dependent steepening of the density spectra around 
SAx marked by the '+'s in the right panel. The vertical dashed line corresponds to the 
outflow wave number /C and all quantities are scaled to the outflow scales. 
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Fig. 5. — Plots showing a slice through log density (top) and log column density (bottom) 
for run HDOI in which turbulence is driven by outflows and isotropic forcing (left panels) 
and for run HDI in which turbulence is driven by isotropic forcing only (right panels). The 
density plots are scaled to p^, and the column density plots to A^o = Po^box- The x and y 
axis are in units of the outflow length scale C 
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Fig. 6. — Histograms of line of sight velocity for runs HDI, HDOI, & HDO weighted by 
density p (upper left), by volume/unweighted (upper right), by line of sight momentum pvz 
(lower left), and by pvl (lower right). Note the high velocity components present in runs 
with outflow forcing (HDOI and HDO) that contribute to the momentum and energy, but 
produce a very small signal in the density weighted emission. 
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Fig. 7. — First three principal eigenimages constructed from PCA analysis of synthetic data 
cubes (left to right) of runs HDI (top), HDOI (middle), and HDO (bottom). All eigenimages 
are normalized to mean emission (column density). Note the lack of large scale coherent 
structures in the first principal eigenimage for runs HDOI and HDO where outflows are 
present, and the large dipole structure present for runs HDOI and HDI where large scale 
forcing is present. 
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Fig. 8. — First three principal eigenvectors constructed from PCA analysis of synthetic data 
cubes (left to right) of runs HDI (top), HDOI (middle), and HDO (bottom). Note the 
broader eigenvectors for runs HDOI and HDI where isotropic forcing is present compared to 
run HDO where only outflow forcing is present. All eigenvectors are normalized to unity. 
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energy produced by large scale isotropic forcing. 



Since outflows drive motions on sub-parsec scales, it seems likely that realistic outflow 
driving will take place in the presence of a externally driven cascade, due to external energy 
inputs into molecular cloud systems from much larger scales. With this motivation, we drove 
outflows into a medium also supplied with isotropcially driven turbulence from larger scales. 
We found that the properties of the turbulence at sub-outflow scales were largely unchanged 
by the external isotropic cascade for typical cloud parameters. In particular, the velocity 
spectra still contain a knee at the outflow scale and a steep slope at sub-outflow scales. We 
also found that the outflows produced a flatter density spectra even at supra-outflow scales. 
We attribute this to the disruption by outflows of larger scale stru ctures produced by the 



large scale driving. We conclude therefore, as in IWang et al.l (|2010[ ). that outflows have the 
potential to play a large role in setting the star formation rate and the initial mass function 
in their neighborhoods. 

Finally, we also find that the density-velocity distributions and to a smaller extent cor- 
relations present in outflow driven turbulence are quite different from those present in an 
isotropic cascade. Many of the techniques for inferring the velocity spectrum of turbulence 
in star forming regions rely on assumptions about density-velocity distributions and corre- 
lations based on simulations of isotropically forced turbulence. To illustrate the resulting 
potential bias, we tested one such model (PCA) against the results from our simulations 
of outflow driven turbulence. We find that indeed PCA underestimates the contribution of 
outflows to the velocity dispersion, though not necessarily the mean scalar momentum. We 
find that PCA is able to infer the driving scale of the outflows for the case of pure outflow 
driven turbulence, but that the presence of even small velocity large scale flows quickly leads 
to an overestimation. We suggest weighting the data cubes by the line of sight velocity as 
one way to make PCA more sensitive to higher velocity flows. This highlights that nature 
of the forcing of turbulence (e.g. isotropic on large scales vs. outflow driven) is of direct 
importance for interpreting observations. 
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Fig. 9. — Plots showing the hne-widths and correlation lengths for the spectroscopically 
resolved principal components. Left panel shows the results for runs HDI, HDOI, and HDO. 
Right panel shows how the line-widths and correlation lengths change when a sinusoidal 
perturbation of varying strength is added to the final frame of run HDO. 
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Fig. 10. — Plot showing how the second principal eigenspectra (top) and eigenimage (bot- 
tom) change when a sinusoidal perturbation of varying strength is added to the final frame 
of run HDO. Plots are scaled as in figures [7] and [HI Note the slight broadening of the eigen- 
spectra and the appearance of a large dipole structure in the eigenimages resulting in a larger 
correlation length {I2) 



